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ALGORITHMS FOR THE DETERMINATION OF THE TIME 
DELAYS OF THE SIGNAL WHEN USING UNEQUAL DETECTORS 

B. L. Novak, 

USSR Academy of Sciences, Institute of Space Research 

In determining the time delays in the reception of /2* 

signals from detectors located in space, it is often advis- 
able to take account of the circumstance that the detectors 
differ both in sensitivity and in time resolution. 

The present work contains a method for determining delays 
which was developed for just such a case. The method is based 
on minimizing the sum of the squares of the discordances of 
the signals recorded by different detectors. The results of 
an investigation of this method are presented by the method 
of maximizing the correlation of the signals. The computed 
algorithms for the evaluation of the accuracy of the solution 
are also presented, once obtained by the means shown in the 
work. 


1. There is a whole series of well-known problems whose /3 
solutions require the determination of the time delays of 
signals detected by several receivers. A related problem is 
the localization in the celestial sphere of the sources of 
gamma-ray bursts [1] . 

In reference [2] the method for determining delays and 
the procedure of evaluating accuracy were developed for the 
case of equivalent detectors. As in [2] , we will deal here 
with detectors having discrete signals (counters) , but, in 


* Numbers in the margin indicate pagination in the foreign text. 
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contrast to [2], we shall take into account that the detectors 
may vary in sensitivity (that is, in the probability of record- 
ing) and in time resolution. 

Let I ( = 1 ,2 ; y = 1 , 2, . . represent the sequence 

of readings of the i^^ detector with a time resolution 
At^^) . We will interpret this sequence as generated by a 
random Poisson process with several variable intensities 
x(i)(.) which are unknown beforehand. It is clear that the 
results produced by the functions assumed above will 

differ for different i in displacement of the argument and 
in factors; that is, there also exist a certain range of t 
and 6 such that all values fulfill the equality 

X(u +T) • (1) 

Frequently the range of t can be limited to a certain region 

! ■ . ■ 

I ITI < To » ^2) 

where tq is a known value. We shall assume the existence of /4 

just Such a situation. 

The problem of finding the time delays therefore reduces 
to determining the values of x satisfying equation (1) by using 
the measurements unknown 3 and fulfilling 

equation (2) . 

Let us examine the possible approaches to the solution to 
the stated problem. We notice that any method must be based 
on a comparison of the measurements and for dif- 

ferent values of the indices, and this reduces to the range 
of those values Z* of the index A at which the sequences 
and correlate, in some sense, best with each other. 


2 



Here the value of the delay which is sought is obviously 
determined exactly at the value At^^) by the expression t* = 

£*At(2) . 


We will not yet take into account the possibility of 
different time resolutions in the detectors, but assume here 
that At(2) = At^^) = At, Mi = M 2 = M, k = (where [•] 

designates the whole integer part of any number) , and m = 

M - 2k. 


Since equation (1) contains the unknown factor 3, it 
is natural to use the method of maximum correlation over the 
range of t. Let the index £ vary from 1 to 2k + 1. Desig- 
nate the coefficient of correlation corresponding to a given 
£ by r£. Then 


rrt 



If max r^ = rg, then the solution to the problem using /5 

the method of maximum correlation is x = £At. 

The requirement of evaluating the accuracy of the solution 
arrived at, however, complicates the use of this method. In 
fact, since r^ is the realization of a random value, in order 
to evaluate the accuracy of the solution one must. know its 
probability value. Analytically obtaining such characteristics 
is a very complicated problem, however. A study of the distrib- 
ution of r^ by means of statistical modeling is also obviously 
unrealistic because of the large (2M) number of parameters to 
be defined. These circumstances force us to turn to other 
methods of locating x. 
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If we try to find the solution to the problem by the 
method of least squares, this reduces to the minimization 
in the given case of the expression 

with respect to Si and: S . Let us first prove that the solution 
to the problem (1) , (2) by the method of least squares is 

equivalent to that of the method of maximum correlation. 

That is, the same value of the index Si represents the maximum 
and minimum of the respective expressions (3) and (4). 


At the maximum of the range of ) the .condition' is 


m 


V 0 ) < 2 ) * -1 r* 

2j riZj "T.® ♦ 


Whence 
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2 0 ) ( 2 ) 
7=1 
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j - « 


Substituting this expression into (4), we arrive at the problem /6 
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The latter is equivalent to the problem 
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V (1) ja) 

■ y--i 


m 


J-i 7=1 

Comparing this expression with (3 ) , we may conclude that the 
value of the index I which represents a minimum in (4) is 

/s 

equal to £. 

The use of the method of least squares makes it possible 
to evaluate the accuracy of the solution in the following 
manner; replace .ji, in (4) with the expression 


ma3C . 
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rn 


2 ^ O) (2) 


7=1 
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r«)» 


2 (2 

and base the subsequent procedure of evaluating accuracy on 
the assumption that 

I ' ■ m 




7-1 

m 

2 -51? 

7=1 


/7 
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We note that it .is also reasonable to represent 6 with 
5-t;" 

the expression - g , ,r . Hereafter we will assume that one 




of these substitutions' has been performed. 


Let us now consider the effect of differences in the time 
resolutions of the detectors on the procedure of determining 
the delays and the evaluation of the accuracy. We assume 
that At^^^ = rAt(2) (5) , where r is an integer. We note 
that the more general case, where r is not an integer,: can be 
reduced to the procedure indicated above. Henceforth Mi 
shall be designated M and Mi-2k as m. We designate 
as ngj for j = 1,2,..., m. 


It is . clear that, in order to satisfy relationship (5) 
and to use the method of least squares, the procedures indi- 
cated above must be based on the minimization of, rather 
than the expression (4) , of the expression 
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2. We now pass to the problem of evaluating the accu- 
racy of the solution. For subsequent analysis it is desirable 
to have the normalized random values which, as in reference 
[2], we obtain by normalizing the sequences andj {^fj] 

by multiplying by and / n (. , respectively. 

Here nj ( .iC t^e/ ^ value of nj and . is the 

averaged value of n^ . 

j 


Using these normalized quantities, we then write 



( 6 ) /8 
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and adopt its use in studying probability characteristics. 


We shall designate random values with capital latin 
letters. Let the parameter of the Poisson value Nj be 
equal to the parameter of be equal to gyj (1) , 


and the parameter of Nv. be equal to y 


: 


‘ p . • 
£>A/j =Py : 


Consequently 


SD^f. 

V 




and 


Whence 


/r*M 
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a-p= = /y^.” ; so -=^ - I 

j/^) ^ ^ 


frij 
✓ft ^ej _ 

l/ A . 


l/^ 


SD 
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= y ?.:- 


and the value x(^)» defined by equation (6), is the realization 
of the random values with a non-central chi-square distribution 
with m degrees of freedom and the parameter of non-centrality: 


■■ ■ . /=»■ 


The value of the parameter can serve as a measure of the 
accuracy of the determination of x. Indeed, we will assume 
the relationship 


( r< (^”‘t l) At 


.yC2). 


is fulfilled for I* satisfying a 


Let us now examine the consequences implied by the use 
of counting devices. For this determination we assume that 


/9 
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such a device is equipped with a second detector and that the 
coefficient of counting is k. 




give an estimate for not for g. Therefore 



just as equation (6), is the realization of the random values 
with a non-central chi-square distribution with m degrees of 
freedom. Thus, the use of counting devices leads merely to 
the appearance of the factor 1/k in the second element of the 
denominator of the coefficient in front of the summation sign 
in the expression for the quantity xC'^)* Henceforth, we will 
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use expression (7) for x(^)» assuming k equivalent units in 
the absence of a counting device. 

Let us return to the problem of evaluating the accuracy /lO 

of the solution. As in reference [2], we will use two methods 
to evaluate the accuracy of the determination of x. The pro- 
cedure of the first method requires an. insignificant replacement 
in the method of the reference mentioned. The interval of 
confidence in is determined by the formula 

h; "=xro-A^+2e“^t24^ 

where ; ic^’J are the upper and lower limits of confidence of 
the non-centrality parameter, is a quantum on the order 

of (2-a)/2 of the normal distribution, and 1-a is the coeffi- 
cient of confidence. 

We now form the system of intervals 
and examine , r which is the union of all those Jg 

w 

for which the condition 

(9) 

is fulfilled. 

The error, Ax, when this method of evaluating accuracy 
is used, is determined by the condition of affiliation of x to 
the expression 

The calculation formulas of the second method for evaluating 
the accuracy of the determination of x are substantially more 
complicated than those in reference [2] . We introduce the 
designations : 
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Using these quantities, the quantity rewritten 

in the form 
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rjO) <l) 

P n 

T 

„fl) (I) 
P n 

/II 

^ V ^ — a 

-{i 





In addition, we introduce the designations: 
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Since p(1-)q( 1) = E, where E is the unit matrix, in the new 
terminology 
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So we obtain for x(^) the expression 


x(e) 


1 + 


ft Af n 


where 
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is the realization of the random normal vector with components /12 

of unit dispersion. Henceforth, elements of the diagonal 

matrix ^ will be designated as p.^ = l,2,...,m, and the 

C 2 ) 

elements of the diagonal matrix Q will be designated = 

1,2, .. . ,M. 


The evaluation of the accuracy in determining t is then 
performed using the statistic 




oT . O o ^ o 

n AgTi - ri Ag n 


/ tn-H . 
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( 10 ) 


where and agg^^^ are diagonal elements of the matrices 

and A^, respectively. Let be a quantum on the order of 
1-a of the normal distribution, d is the error of the second 
order, and Jg, as in the first method of evaluating accuracy, 
is the interval 


If 


( 11 ) 


11 



then the corresponding interval Jg will be included in the 
union , and the error At is determined by the condi- 

tion of affiliation of t to the expression 

3. Since the dimensions of the matrices and in 

expression (10) - 2M - are very large mambers, operating with 
these matrices immediately is difficult. In this connection 
we can obtain simple expressions describing the dependence of 
the statistic on the sequential measurements | j 

We designate as and itiax f t, t] as Izf so that 

£2 ^ £ 1 . Clearly, 



If the non-zero blocs of this matrix are designated 


Af-Af = 


A„ 

- A-f 

£ 

£ 





A'"'' 



is rewritten in the form 


And if we introduce the designation V for , and W for 

^ then the quadratic form from 

expression (10) is obviously written in the form + ww"'^ . 

The formulas for the calculation of this quadratic form depend 
on the relationships among £ 1 , £ 2 , M, m and r. 


Let ^0 = 
number , and 
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where [•] is the integral part of the 
£ 1 ~ £ 0 • 


1) If the relationship £ 1 + rm - 1 ^ £2 is fulfilled, 

then 
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2) If the relationship il i + rm - 1 < £2 is fulfilled. 


then 
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Let us pass to the calculation of Sfi (A^ - 4/)® :: from expres- 
sion (10) . As is well known, any quadratic matric B of dimen- 
sion n with elements 6.. fulfills the condition 

V 

It follows from this equation that 

. ^ fn*M. M-.' 
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Here 6aj^j are the elements of the matrix - Aj, a^j are 
elements of the matrix A^^^, and are elements of the 


matrix A 
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We obtain from the last equation the following: 
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Finally, for the final factor in brackets in expression (IQ) 
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the correct formulas are 
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Thus, the algorithm of the accuracy of the second methods 
consists of calculating the statistic (10) for all values of 
I, using formulas (12), (13), (16) and (18) for 5 , 1 + rm ^ £2 
and (14), (15), (17) and (19) for £ 1 + rm < £ 2 , while checking 
that condition (11) is fulfilled. 
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